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Abstract. By means of numerical simulations we investigate the configurational properties of densely and 
fully packed configurations of loops in the negative-weight percolation (NWP) model. In the presented 
study we consider 2d square, 2d honeycomb, 3d simple cubic and 4d hypercubic lattice graphs, where edge 
weights are drawn from a Gaussian distribution. For a given realization of the disorder we then compute 
a configuration of loops, such that the configurational energy, given by the sum of all individual loop 
weights, is minimized. For this purpose, we employ a mapping of the NWP model to the "minimum-weight 
perfect matching problem" that can be solved exactly by using sophisticated polynomial-time matching 
algorithms. We characterize the loops via observables similar to those used in percolation studies and 
perform finite-size scaling analyses, up to side length L — 256 in 2d, L = 48 in 3d and L = 20 in 4d (for 
which we study only some observables), in order to estimate geometric exponents that characterize the 
configurations of densely and fully packed loops. One major result is that the loops behave like uncorrelated 
random walks from dimension d = 3 on, in contrast to the previously studied behavior at the percolation 
threshold, where random-walk behavior is obtained for d > 6. 



1 Introduction 

■ The statistical properties of lattice-path models on graphs, 
I equipped with quenched disorder, have experienced much 
, attention during the last decades. They have proven to 

■ be useful in order to describe line-like quantities as, e.g., 
, linear polymers in disordered media [1 2 3 4j . vortices in 

■ high Tc superconductivity [516] , cosmic strings in the early 

■ universe |7I8I9| , and domain- wall excitations in disordered 
. media such as 2d spin glasses [1011 and the 2d solid-on- 
solid model [12 ■ The precise computation of these paths 
can often be formulated in terms of a combinatorial opti- 
mization problem and hence might allow for the applica- 
tion of exact optimization algorithms developed in com- 
puter science |13I14I15] . So as to analyze the statistical 
properties of these lattice path models, geometric observ- 
ables and scaling concepts similar to those used in per- 

■ eolation theory |16I17| or other "string" -bearing models 
' |18ll9j are often applicable. 

The topic of the presented article is the negative-weight 
percolation (NWP) problem |20I21I22] . wherein one con- 
siders a regular lattice graph with periodic boundary con- 
ditions (BCs) and where adjacent sites are joined by undi- 
rected edges. Weights are assigned to the edges, represent- 
ing quenched random variables drawn from a distribution 
that allows for edge weights of either sign. For a given real- 
ization of the disorder, one then computes a configuration 
of loops, i.e. closed paths on the lattice graph, such that 
the total sum of the weights assigned to the edges that 
build up the loops attains an exact minimum. Note that 



the application of exact algorithms in contrast to standard 
sampling approaches like Monte Carlo simulations avoids, 
e.g., equilibration problems. Also, since the algorithms run 
in polynomial time, large instances can be solved. As an 
additional optimization constraint we impose the condi- 
tion that the loops are not allowed to intersect. Conse- 
quently, since a loop does neither intersect with itself nor 
with other loops in its neighborhood, it exhibits an "ex- 
cluded volume" quite similar to usual self avoiding walks 
(SAWs) [17]. The problem of finding these loops can be 
cast into a minimum-weight path (MWP) problem, out- 
lined below in sect. [5] in more detail. 

In previous studies [20121122] , the details of the weight 
distribution were further controlled by a tunable disor- 
der parameter. As a pivotal observation it was found 
that, as a function of the disorder parameter, the NWP 
model features a disorder driven, geometric phase transi- 
tion. This transition leads from a phase characterized by 
only "small" loops to a phase that also features "large" 
loops that span the entire lattice along at least one direc- 
tion. Regarding these two phases and in the limit of large 
system sizes, there is a particular value of the disorder pa- 
rameter at which percolating (i.e. system spanning) loops 
appear for the first time [3D]. Previously, we have inves- 
tigated the NWP phenomenon for 2d lattice graphs [5D] 
using finite-size scaling (FSS) analyses, where we char- 
acterized the underlying transition by means of a whole 
set of critical exponents. Considering different disorder 
distributions and lattice geometries, the exponents were 
found to be universal in 2d and clearly distinct from those 
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Fig. 1. Illustration of densely packed loop (DPL) config- 
urations, (a) Sample configuration of DPLs for a square 
lattice of side length L = 64 and fully periodic bound- 
ary conditions (BCs), and (b) configuration of DPLs for 
a simple cubic lattice of side length L = 32 and fully peri- 
odic BCs. For a more clear presentation, small loops with 
length ^< 10 and a very large loop with length i — 26 770 
were omitted from the latter figure. 



describing other percolation phenomena. In a subsequent 
study we investigated the effect of dilution on the critical 
properties of the 2d NWP phenomenon 21 . Therefore we 
performed FSS analyses to probe critical points along the 
critical line in the disorder-dilution plane that separates 
domains that allow/disallow system spanning loops. One 
conclusion of that study was that bond dilution changes 
the universality class of the NWP problem. Further we 
found that, for bond-diluted lattices prepared at the per- 
colation threshold of 2d random percolation and for purely 
Gaussian distributed edge-weights, the geometric proper- 
ties of the system spanning loops compare well to those 
of ordinary self-avoiding walks. Lately, we performed fur- 
ther simulations for the NWP model on hypercubic lat- 
tice graphs in dimensions d = 2 through 7 "221 , where we 
found evidence for an upper critical dimension du — 6 of 
the NWP phenomenon. To resume, up to now we have fo- 
cused on the critical properties of the NWP model in the 
vicinity of the critical point at which percolating loops 
first appear in the limit of large system sizes. As we expe- 
rienced, at this critical point the loops are rather isolated 
and well separated from each other, i.e. they resemble a 
dilute gas of loops. However, as the disorder on the lattice 
increases, the loops get more dense and it is reasonable to 
ask for the statistical properties of these densely packed 
configurations of loops. 

Here, we investigate the NWP model on 2d square 
and honeycomb, Sd simple cubic and Ad hypercubic lat- 
tice graphs for a purely Gaussian distribution of the edge 
weights. Therein, our aim is to characterize the minimum- 
weight configurations of densely packed loops (DPLs) , see 
Fig. [U regarding their statistical properties. Further, we 
present a modification of the NWP algorithm that allows 
for the computation of configurations of fully packed loops 
(FPLs). Regarding a configuration of FPLs, each site on 
the lattice is visited by a loop. We find that the scaling 
properties of DPLs and FPLs for the 2d systems agree 
within errorbars and compare well to those obtained for 
FPLs related to other discrete interface models [531. Fur- 



ther, the geometric properties of the loops for Sd systems 
are similar to those of completely uncorrelated random 
walks, indicating that the upper critical dimension shifts 
to d^^^ = 3 for DPLs. 

The remainder of the presented article is organized as 
follows. In section[21 we introduce the model in more detail 
and we outline the algorithm used to compute the mini- 
mum weight configurations of loops (i.e. DPLs and FPLs). 
In section [31 we present the results of our numerical sim- 
ulations and in section |4] we conclude with a summary. 



2 Model and Algorithm 

In the remainder of this article we consider regular 2d 
square/honeycomb, 3d simple cubic and 4d hypercubic 
lattice graphs G={V, E) with side length L and fully pe- 
riodic boundary conditions (BCs). The considered graphs 
have N=\V\=L'^ sites ieV and a number of \E\=zN/2 
undirected edges {i,j}€E that join adjacent sites i,j£V. 
Above, z signifies the coordination number of the lattice 
geometry, where z = 4, 3, 6 and 8 for square, honey- 
comb, simple cubic and Ad hypercubic lattice graphs, re- 
spectively. We further assign a weight uJij to each edge 
contained in E. These weights represent quenched ran- 
dom variables that introduce disorder to the lattice. In 
the presented article we consider independent identically 
distributed weights drawn from a Gaussian disorder dis- 
tribution with mean zero and unit width, i.e. 



P{uj) = exp (-w72)/V27r. 



(1) 



Note that this distribution explicitly allows for loops £ 
with a negative total weight loc = "^u aec ^y • ^^ ^^ ^^ 
support intuition: Eq. [T] corresponds to the limit p ~ 1 
of the disorder parameter for the edge-weight distribu- 
tion considered in previous studies of the NWP model, 
see e.g. Ref. [22]. Accordingly, we here consider the NWP 
phenomenon high up in the percolating phase where the 
minimal weight configurations of loops are densely packed. 
The NWP problem then reads as follows: Given G to- 
gether with a realization of the disorder, determine a set 
C of loops such that the configurational energy, defined 
as the sum of all the loop- weights £ = X^rec '^-C' i^ niin- 
imized. As further optimization constraint, the loops are 
not allowed to intersect. Clearly, the configurational en- 
ergy £ is the quantity subject to optimization and the 
result of the optimization procedure is a set of loops C, 
obtained using an appropriate transformation of the orig- 
inal graph as detailed in [21]. For the transformed graphs, 
minimum-weight perfect matchings (MWPMs) [25126127] 
are calculated, that serve to identify the loops for a given 
realization of the disorder. This procedure allows for an 
efficient implementation |23 of the simulation algorithms. 
Until further notice, the subsequent description applies to 
DPLs only. In order to obtain FPLs, the transformation 
is slightly different (see discussion below). Here, we give a 
brief description of the algorithmic procedure that yields 
a minimum-weight set of loops for a given realization of 
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the disorder. Fig. ^ illustrates the 3 basic steps, detailed 
below: 

(1) each edge, joining adjacent sites on the original 
graph G, is replaced by a path of 3 edges. Therefore, 2 
"additional" sites have to be introduced for each edge in 
E. Therein, one of the two edges connecting an additional 
site to an original site gets the same weight as the cor- 
responding edge in G. The remaining two edges get zero 
weight. The original sites i £ V are then "duplicated", 
i.e. i — >■ ?i,i2, along with all their incident edges and the 
corresponding weights. For each of these pairs of dupli- 
cated sites, one additional edge {11,12} with zero weight 
is added that connects the two sites «i and Z2. The result- 
ing auxihary graph Ga = {Va,Ea) is shown in Fig. [2Ib), 
where additional sites appear as squares and duplicated 
sites as circles. Fig. [2jb) also illustrates the weight assign- 
ment on the transformed graph Ga- Note that while the 
original graph (Fig. [H^a)) is symmetric, the transformed 
graph (Fig. [2l^b)) is not. This is due to the details of the 
mapping procedure and the particular weight assignment 
we have chosen. A more extensive description of the map- 
ping can be found in 11 . 

(2) a MWPM on the auxiliary graph is determined 
via exact combinatorial optimization algorithms |29| . A 
MWPM is a minimum- weighted subset M of Ea, such 
that each site contained in Va is met by precisely one 
edge in M. This is illustrated in Fig.j^Jc), where the solid 
edges represent M for the given weight assignment. The 
dashed edges are not matched. Due to construction, the 
auxiliary graph consists of an even number of sites and the 
transformation procedure described in step (1) guarantees 
that a perfect matching exists. Note that a MWPM can be 
computed in polynomial time as a function of the number 
of sites, hence large systems with several thousands of sites 
are feasible. 

(3) finally it is possible to find a relation between the 
matched edges M on Ga and a configuration of negative- 
weighted loops C on G by tracing back the steps of the 
transformation (1). As regards this, note that each edge 
contained in M that connects an additional site (square) 
to a duplicated site (circle) corresponds to an edge on 
G that is part of a loop, see Fig. [D^d). Note that, by 
construction of the auxiliary graph, for each site ii or 
12 matched in this way, the corresponding twin site 12/?! 
must be matched to an additional site as well. This guar- 
antees that wherever a path enters a site of the original 
graph, the paths also leaves the site, corresponding to the 
defining condition of loops. All the edges in M that con- 
nect like sites (i.e. duplicated-duplicated, or additional- 
additional) carry zero weight and do not contribute to a 
loop on G. Once all loop segments are found, a depth-first 
search |24I26| can be used to identify the loop set C and 
to determine the geometric properties of the individual 
loops. Here, the exemplary weight assignment illustrated 
in Fig. [Sfa) yields 2 loops, i.e. C = {Ci,C2}, with weights 



(a) 



W£i = 0JC2 = -4 and lengths £1 = J2{i,j}eCi 1 = 8, 4=4. 
Hence, the configurational energy reads £ = —8. 

The result of the calculation is a collection C of loops 
such that the total loop weight, and consequently the con- 
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Fig. 2. Illustration of the algorithmic procedure: (a) orig- 
inal lattice G with edge weights. For clarity, a bimodal 
distribution that yields edge- weights ±1 is considered. 
Further, only negative edge- weights are shown. Unlabeled 
edges have weight -1-1. (b) auxiliary graph Ga with proper 
weight assignment. Black edges carry the same weight as 
the respective edge in the original graph and gray edges 
carry zero weight, (c) minimum- weight perfect match- 
ing (MWPM) M: bold edges are matched and dashed 
edges are unmatched, and (d) loop configuration (bold 
edges) that corresponds to the MWPM depicted in (c). 
(e) MWPM for the the modified mapping that yields fully 
packed configurations of loops, (f) fully packed loops (bold 
edges) that correspond to the MWPM shown in (e). 



figurational energy £, is minimized. Hence, one obtains a 
global collective optimum of the system. Obviously, all 
loops that contribute to C possess a negative weight. Also 
note that the choice of the weight assignment in step (1) is 
not unique, i.e. there are different possibilities to choose 
a weight assignment that all result in equivalent sets of 
matched edges on the transformed lattice, corresponding 
to the minimum- weight collection of loops on the original 
lattice. Some of these weight assignments result in a more 
symmetric transformed graph, see e.g. |24) . However, this 
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is only a technical issue that does not affect the resulting 
loop configuration. Finally, for an illustrational purpose, a 
small 2d lattice graph with free BCs was chosen intention- 
ally. The algorithmic procedure extends to higher dimen- 
sions and fully periodic BCs in a straight-forward manner. 

Note that, by means of the algorithmic procedure out- 
lined above, its also possible to compute configurations 
of fully packed loops. To accomplish this, only one detail 
related to step (1) has to be altered. I.e., original sites 
i GV are duplicated (i —>■ 11,12) to yield 11,12 G Va, but 
no additional edge {^1,12} is added to Ea- By construc- 
tion, a MWPM, as computed in step (2), will not contain 
an edge that connects a tuple of duplicated sites (since 
there are no such edges present). Hence, in the resulting 
MWPM, each duplicated site will be matched to an addi- 
tional site, see Fig.j^fe). Applying step (3), one then finds 
that a each site on the original graph is an endpoint of ex- 
actly two loop segments, see Fig. ^f), and the result is a 
fully packed configuration of loops. Applying this modified 
mapping to the exemplary weight assignment illustrated 
in Fig. Ha) yields Cfpl = {Ci,£2} (see Fig. [^f)), with 
weights ojCi = —6, ujc2 — and lengths £1 = 12, £2 = 4. 
Here, the configurational energy reads £ = —6. 

Apparently, there are qualitative differences between 
configurations of DPLs and FPLs: Regarding the DPLs, 
not all lattice sites are visited by a loop (see Fig. da)), 
and, due to the "energy minimization principle" of the 
optimization procedure, the weight of an individual loop 
is strictly smaller than zero. In contrast, FPLs require 
each lattice site to be visited by a loop. This constitutes 
a "hard" constraint for the optimization procedure and 
as a result, an individual loop weight might take a value 
larger than or equal to zero (as for the assignment of ±1 
edgeweights in Fig. [2Ia)). 

Subsequently, we will use the procedure outlined above 
so as to investigate the NWP model on 2d square and hon- 
eycomb and 3d simple cubic lattice graphs with a Gaussian 
distribution of edge weights. 



3 Results 

In order to characterize the statistical properties of config- 
urations of DPLs we performed simulations for 2d square 
{2d-sq) systems with a linear extension of up to L = 256 
lattice sites, and 3d simple cubic (3d-sc) systems with side 
length up to L = 48 sites. If not stated explicitly, the results 
presented in the remainder of this section follow from an 
analysis of the full ensemble of loops obtained for a fixed 
system size, i.e. L = 256 (48) in 2d (3d). For both setups 
we considered a number of n= 12800 independent realiza- 
tions of the disorder to compute averages. Moreover, we 
conducted simulations for a 2d honeycomb (2d-hy) DPL 
and 2d-sq FPL version of the loop model for lattices up 
to L = 128, where we considered n = 9600 samples for 
the largest systems. We performed additional simulations 
for Ad hypercubic (4d-hc) lattice graphs with L = 20 and 
n — 9600 in order to study selected observables for an even 
higher dimension. 



Next, we briefly introduce the basic observables we 
used to characterize the individual loops obtained from 
the simulations. As introduced earlier, two fundamental 
quantities related to a given loop are its weight ujc and 
length £. We further determine the linear extensions i?i, 
i — 1 . . .d, of each loop by projecting it onto the inde- 
pendent lattice axes (for the hexagonal lattice geometry 
we therefore use the topologically equivalent brickwall lat- 
tice). The largest of those values, i.e. 



R- 



niax (Ri), 

1=1. ..d 



(2) 



is referred to as the spanning length of the loop and the 
smallest, i.e. 



r— niin (Ri), 

i=l...d 



(3) 



is called its roughness. Moreover, so as to characterize the 
full perimeter of an individual loop on a coarse grained 
scale, we can define the size 



Rs = 2_^ Ri 



(4) 



j=i 



i.e. the length of the loop if all small scale irregularities 
were flattened [7]. Note that the linear extensions of a 
given loop are measured in units of lattice sites. E.g., the 
loop £1 contained in the example illustrated in Fig. [IJd) 
has i?i = i?2 = 3 and thus i?s = 6. Albeit this definition of 
a loop size is somewhat odd for small loops in a 3d setup 
(where an elementary loop of length £ = 4 has size i?s=5, 
as noted in Ref. [7]), it nevertheless gives a reasonable 
deflnition of a coarse grained loop size in the limit of large 
£. Anyway, note that there are some difficulties related 
to the "extremal" loops on the lattice. This means, if we 
consider the full ensemble of loops for a system of size L, 
there are very small loops {R <^ L) which are affected by 
the finite lattice spacing and also very large loops (i?«L) 
that are affected by the finite size of the considered system. 
So as to obtain most reliable estimates for the observables 
that characterize the configurations of loops, and if not 
stated explicitly, these extremal loops are discarded from 
the analysis. Subsequently, we will adopt the notation (O) 
to signify the disorder averaged value of observable O. 
Whenever a more clear presentation requires, we will write 
{0)x to signify the disorder averaged value of O with 
respect to the quantity X . 

In the following we will present the results on the geo- 
metric properties of the DPL/FPL configurations in sub- 
section [5TTl the results on the energetic properties are sum- 
marized in subsection [ 



3.1 Geometric properties 

As pointed out above, we here report our results regard- 
ing the geometric properties of the loop configurations ob- 
tained for the 2d (3d, 4d) systems mention in the beginning 
of the presented section. 
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Table 1. "Packing" properties, scaling exponents and parameters that characterize DPLs/FPLs on the considered 
lattice graphs. From left to right: System setup [sq=square, hy^honeycomb, sc^simple cubic, hc=hypercubic; DPLs 
(FPLs)=densely (fully) packed loops], linear extension L of the lattice for which the listed estimates were obtained, 
number n of samples considered, approximate number Moops of loops collected during the simulations, probability 
pcov that a lattice site is visited by a loop, fraction foe of loop segments contained in system spanning loops and 
"loop-shape" parameter a/{2(P). Fractal dimension df, correction to scaling exponent uj and exponent r that describe 
the scaling of the non-spanning loops. Scaling dimension d'^ for the largest loops found for each realization of the 
disorder and asymptotic estimate e for the normalized configurational weight. 







L 


n 


-i^loops 


pcov 


/ex. 


a/{2S) 


df 


LO 


r 


d'f 


(^> 


2d-sq 


DPLs 


256 


12 800 


25 X 10^ 


0.825 


0.327 


0.97(3) 


1.756(8) 


0.86(8) 


2.14(2) 


1.752(1) 


0.596 


2d-sq 


FPLs 


128 


9 600 


6x 10® 


1.000 


0.371 


0.98(1) 


1.75(2) 


0.9(2) 


2.16(1) 


1.751(1) 


0.493 


2d-hy 


DPLs 


128 


9 600 


2 X 10'' 


0.750 


0.429 


1.04(4) 


1.75(2) 


1.0(1) 


2.14(1) 


1.750(1) 


0.412 


3d-sc 


DPLs 


48 


12 800 


13 X 10** 


0.901 


0.907 


0.96(1) 


2.01(4) 


0.66(11) 


2.49(6) 


3.00(1) 


0.844 


4d-hc 


DPLs 


20 


9 600 


6x lO** 


0.938 


0.970 


0.99(3) 


2.00(5) 


- 


- 


- 


1.013 



Global "packing" properties: As pointed out above, for 
the configurations of DPLs considered here, not all lattice 
sites are visited by a loop. We can quantify this by comput- 
ing the probability pcov that some site on the lattice lies on 
the perimeter of a loop, see Tab. [1] Via additional simula- 
tions for other system sizes (not shown here), we verified 
that the value of pcov is basically independent of the lat- 
tice sizes. From the accumulated data we also obtained 
the probability poo that some site on the lattice lies on 
the perimeter of a large, i.e. system spanning, loop. From 
this we conclude that a fraction foo^ Poo/ Pcov (see Tab.[T]) 
of all loop segments for the 2d-sq, 2d-hy, Sd-sc and Ad-hc 
setup is contained in large loops. Further, note that the 
estimate /oo ~ 0.371 for the 2d-sq FPL version of the loop 
model (where pcov = Ij exactly) is somewhat larger than 
the corresponding value foo « 0.327 for the DPL configu- 
rations. As evident by comparing the estimates pcov and 
/oo obtained for DPLs on 2d square and honeycomb lattice 
graphs, the observables above characterize the particular 
lattice geometries chosen and hence are not universal. 



Scaling of the loop shape: Next, we want to study the 
scaling behavior of the average loop shape. As in Refs. 
[7122] we therefore monitor the volume to surface ratio 
Vb/Sb of the smallest box that fits the individual loops 
as a function of the coarse-grained loop size Rs, where 
Vb = nf=iRz and S'b = 2 X Y.Li Vb/Ri- For hypercubic 
volumes with identical values Ri — . . . ^ Rd, one would 
expect to find Vb/'S'b = (2(i^)"^-Rs. Performing fits to the 
form {Vb / Sb) = aRf , discarding very small and large non- 
percolating loops, we here yield the results V' — 1.002(5) 
(1.001(1)) and a = 0.121(3) (0.0535(3)) that describe 2d- 
sq (3(i-sc) configurations of DPLs. Further, for the 2d-sq 
FPL data we obtain tp = 0.999(2) and a = 0.123(1). For 
the 2d-hy setup our data is consistent with the precise 
scaling expression {Vb/Sb) = 0.130{5){Rs ~ 0.4)ioo6(8)^ 
that includes corrections to scaling in a very basic form. 
The same holds also for the 4d-hc setup, where (Vb/S'b) = 
0.031(l)(i?s-0.5)°-994(s). Hence, we find values of a/(2d2) 
(see Tab.d]) reasonably close to 1 in order to conclude that. 



in a statistical sense, the loops are not oblate but possess 
a rather spherical shape. 



Scaling of the loop length - "small" loops: The geomet- 
ric properties of the loops can further be characterized by 
the fractal dimension d^. It can be defined from the scal- 
ing of the average loop length {€) as a function of the loop 
size according to df = limfl^_>oo log((^))/log(i?s)- Lower 
and upper bounds on the above estimate of dj can be 
obtained in a systematic way [30]: Based on the expres- 
sion Rs = 'Ylii=\ ^i for the loop size, it is intuitive that 
Rj = d ■ r typically underestimates the value of Rs (bear 
in mind that r ~ mini=i...rf(i?i)). Consequently, the scal- 
ing assumption {£)r ^ r f yields an estimate d7 < df. In 

the same sense, Rf = d ■ R overestimates Rg and gives 
rise to a scaling exponent dt > df. However, in the limit 
of large loops and since the loops have rather spherical 
shape (which implies r ^ R), we should be able to verify 
dj = df = dt . For the 2d and 3d systems considered, we 
found our data best fit by a scaling function that includes 
corrections to scaling in the form {£) ^ X'^f {1 + cX~'^), 
where X optionally stands for r, Rg or R. In this regard 
we yield the estimates listed in Tab. [5] and fit-parameters 
c of the order of unity (Note that for the 4(i hypercube 
with L = 20, since the range of values for the scaling pa- 
rameters r and R is too small, our data did not allow for 
a decent analysis of upper and lower bounds). The expo- 
nents for the choice X — Rs are listed in Tab. [T] and the 
data for the 2d-sq (3d-sc) DPLs is shown in Fig. [31 When 

comparing the resulting estimates of dJ and df, keep 
in mind that the presence of corrections to scaling hinders 
the straight forward fit to a pure power law function. As 
an alternative to the scaling function considered above, 
one might choose a more basic expression that accounts 
for corrections to scaling and features less free parame- 
ters. In this regard, for the 2d-sq FPL setup we obtained 
{i) ^ {Rs + 1.5(4))-'^-^'*^*^^\ in agreement with the estimate 
listed in Tab. [TJ From the simulations for the 4d systems 
and in the limit of large loops, we found a good fit of 
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2(/-DPLs: L=256 
3(/-DPLs: L=48 
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Fig. 3. Scaling of the average loop length {£) of the non- 
spanning DPLs as a function of the coarse grained loop 
size Rs for 2d square and Sd simple cubic lattice graphs 
(the Sd data is shifted upwards by a factor of 3). Further, 
the upmost data curve indicates the scaling ~ i?^ of Sd 
loops with spanning index j = l, i.e. those loops that span 
the system along exactly one direction (the respective data 
points are shifted upwards by a factor of 12). The inset 
illustrates the sequence of effective scaling exponents d^I^ 
as function of the inverse loop size for the DPL setup on 
the 2d square lattice. 



the data to the form {£) ^ Rs^ , resulting in the estimate 
d/ = 2.00(5). 

As further check of the results for the 2d-sq (3d-sc) 
DPLs obtained above, we computed the effective (lo- 
cal) slopes that describe the scaling of {£) within inter- 
vals of, say, m consecutive data points Rg according to 
log((£)) = dj^log(i?s) -I- c. By sliding this "scaling win- 
dow" over the whole data set, one obtains a sequence of 



Table 2. Extended analysis of scaling exponents and pa- 
rameters that characterize DPLs/FPLs. From left to right: 
lattice setup [sq=square, hy=honeycomb, sc=simple cu- 
bic; DPLs (FPLs) = densely (fully) packed loops], lower 

and upper bounds d7 on the fractal dimension and cor- 
responding correction to scaling exponents w"*"' that de- 
scribe the scaling of the non-spanning loops as explained 
in the text. 



d- 



U] 



4 



2d-sq DPLs 

2d-sq FPLs 

2d-hy DPLs 

3d-sc DPLs 



1.749(6) 
1.73(2) 
1.73(4) 
2.0(1) 



0.8(3) 
1.1(3) 
1.0(3) 
0.5(2) 



1.751(5) 0.92(5) 



1.75(1, 
1.75(2) 
2.0(3) 



0.7(1) 
1.1(1) 
0.3(4) 
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Fig. 4. Scaling behavior of the probability density func- 
tion Pl {(-) of the loop length I, considering the largest loop 
found for each realization of the disorder for different sys- 
tem sizes L. After rescaling, data curves for the different 
2d setups [DPLs (FPLs) = densely (fully) packed loops; 
sq=square, hy=honeycomb] fall onto one universal scaling 
function. Further, the contribution of the loops with span- 
ning index j = 0, 1 , 2 (see text) to the L = 256 (2d DPL-sq) 
data curve is illustrated. 



effective exponents d%^{Rs). The asymptotic value d°l^ is 
then extrapolated by means of a straight line fit to the 
plot of d'i^(i?s) as a function of the inverse loop size. For 
the 2d-sq (3d-sc) DPL data, setting to = 20 (10), this pro- 
cedure yields the estimates df = 1.749(6) (1.98(3)), that 
agree with those listed in Tab. [T] within errorbars. The 
results for the DPLs on the 2d square lattice are further 
shown in the inset of Fig. [31 

Referring to the numerical values of d/ found here, and 
from a point of view of statistical physics, the loops appear 
to be self similar fractals. Albeit the loops are self avoiding 
by construction, our results indicate that for d > 3 they ex- 
hibit the statistical properties of completely uncorrelated 
random curves with df = 2. Hence, this suggests that the 
upper critical dimension of the NWP problem in the limit 
of DPLs is d!?PL = 3. 



Scaling of the loop length - "large" loops: So far, the 
data analyses referred to the subset of non-spanning loops 
found for the simulated disorder samples. However, as dis- 
cussed above, about 91% of the over all loop length for 
the 3d-sc systems is contained in very large, i.e. system 
spanning, loops (see Tab.[T]). For this setup we further ob- 
serve that, in comparison to the small loops investigated 
above, very large system spanning loops exhibit different 
scaling properties. These loops are a direct consequence 
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of the periodic boundary conditions, and, although they 
are strongly affected by the finite size of the "simulation 
box" , we might nevertheless characterize their configura- 
tional properties. 

To perform a systematic analysis of these very large 
loops, we here introduce the spanning index j € {0 . . .d} 
of an individual loop. As regards this, a value of, say, j = 2 
indicates that the respective loop spans the lattice along 
exactly 2 of the d independent lattice directions (non- 
spanning loops have j — 0). Note that for a single configu- 
ration of the disorder, the number of loops having j > 1 is 
rather small. E.g., simulating n = 12 800 instances for a 3d 
simple cubic system with L = A8 we count ^loops"^ 149, 
3 260, 31 106 for j — 1, 2, 3, respectively (for comparison 
-^ioops ~ 1^ ^^^ 016). Hence, the results below should be 
recognized as "guiding values" that are to be encountered 
with some caution. 

Now, filtering the 3d data for loops with spanning in- 
dex j = 1, we observe a power law scaling of the form 

(£)~i?s' , see Fig.ini Considering a loop with spanning ,j^, o 

index j = 2, one might argue that the roughness r is a more 
adequate scaling variable. Note that for loops with j = 2 
on a lattice of side length L = 48 one has r = i?^ — 2 • 48. As 
regards this, upon analysis we find a best fit of the data to 
the form {£) ^ {Rs - 91(8))^'^'^^^. However, as mentioned 
above, the number of loops having j = 2 is rather small. 
Hence the statisics for these loops is rather poor and the 
above estimate presents only an approximate result. Note 
that no such scaling analysis is possible for a spanning in- 
dex j = d, since the corresponding loops all have Rs=d- L. 

Regarding the 4d data we yield iVj^^p^ = 6625020, 2873, 
1 939, 2 021 and 30 961 for j = 0, . . . , 4, respectively. Due 
to the small number of loops having j = 1, 2, 3, no decent 
scaling analysis was possible. 

We also performed FSS analyses according to {I) ^ 
L f , where, for different system sizes L, we consider only 
the largest loop found for each realization of the disorder. 
In this regard we were able to verify the "small loop" 
scaling behavior for the 2d systems, i.e. d'f ~ 1.75, see 
Tab. [T] A similar scaling analysis for the 3d-sc data yields 
(^) ^L'^"°'-^^ (not shown), where all loops considered (i.e. 
the largest loop for each realization of the disorder) turned 
out to have j = 3. In Ad we performed simulations for 
L = 20 only. Hence no such scaling analysis is carried out 
for the M setup. 

The probability density function (pdf) Pl{^) of the 
length t for the largest loop for a given realization of the 
disorder is shown in Fig. 01 As evident from the figure, 
the data for all the 2d setups considered fall onto one 
master curve, if rescaled according to Pl {tj = {St)2^po [{£— 
{^) l)I {5()l]- Therein, pq[-] signifies a scaling function and 
{5(.)l = [{(■'^)l — (^)i]^^^ refers to the root mean square 
deviation related to the loop length £. This data collapse 
highlights that the different geometries and types of loops 
exhibit an universal behavior. We further observed that 
each data curve is composed of three parts that relate to 
loops with different values of j. For the 2d DPLs on the 
L = 256 square lattice. Fig. 2] also shows the contribution 
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Fig. 5. Results obtained for the number densities n{-) 
for loops of length (. and roughness r. The main plots 
show the number density n(£) of loops with length £ (using 
logartihmic binning) and the inset illustrates the scaling 
of the number density n{r) of loops with roughness r. (a) 
Data obtained for DPLs on the 2d square lattice, (b) data 
for DPLs on the 3d simple cubic lattice. 



stemming from loops with spanning indices j = 0, 1 and 2. 
Loops with j = contribute to the shortest loops, j = 2 
the largest loops and j ^ 1/ j ~ 2 contribute to the bulk 
of itermediate size loops. 



Scaling of the loop area in 2d: For the 2d setups we can 
further associate a cluster with each loop. Such a cluster 
simply consists of the elementary plaquettes enclosed by 
a loop (for a square lattice, an elementary plaquette con- 
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sists of a unit square bordered by 4 lattice edges) . Conse- 
quently, the loop forms the perimeter of its associated clus- 
ter. The basic observable related to a cluster is its area (or 
similarly its volume) Ac , measured by the number of ele- 
mentary plaquettes that comprise the cluster. To support 
intuition, the exemplary loop configuration C = {£i,C2} 
shown in Fig. [2jd) features two clusters with perimeter 
ii = 8, ^2 = 4 and area Aci = 3, Ac2 = 1. 

As introduced earlier, the spanning length i? is a typ- 
ical lengthscale that characterizes an individual loop, and 
henceforth characterizes also the associated cluster. The 
area of the clusters is expected to scale as Ac ~ -R''" , where 
dy signifies the fractal dimension of the cluster area (or 
volume). For 2d-sq DPLs we find (Ac) -^ R^'^^'^^*^ (2d-hy 
DPLs: d„ = 2.01(1); 2d-sq FPLs: d^, = 1.99(1)), indicating 
that the clusters are compact with dy = d. Regarding the 
fractal scaling of the cluster perimeter {l^W^f), we might 
further write Ac ^^ i'^^'^^' , see Refs. [31132] . More precise, 
to account for corrections to scaling we consider a scaling 
function {Ac)^Ao{£ + A£Y and obtain 





^0 


/It 


K 


DPLs - sq : 


0.378(2) 


-1.58(6) 


1.137(1) 


FPLs - sq : 


0.323(2) 


-0.71(4) 


1.135(1) 


DPLs - hy : 


0.397(9) 


-2.2(2) 


1.135(4) 



Here, since df —d/n, the scaling exponent k provides an 
independent way to estimate the fractal dimension of the 
loops. In any case, the numerical value of d/ k is in good 
agreement with the previous estimate of df for the respec- 
tive setup listed in Tab. [TJ 



Number densities: First, we consider the scaling behav- 
ior of the number density n{Rs) of small loops with size 
Rs- Upon analysis of the 2d and 3d data, we obtain 
scale invariant expressions of the form n{Rs) '^ R7 
(not shown). Argumenting as above, it should not mat- 
ter whether we consider the number density of loops with 
given spanning length or roughness instead of n{Rs). In 
this regard, considering only non-spanning loops, we find 
n{r) ~ r-3™(3) (2d-sq DPLs) and n{r) - r-3-95(6) (3^. 
sc DPLs), see inset of Figs. [51^ a), (b), where in either case 
the algebraic decay is governed by an exponent in good 
agreement with {d+1). Note that, if the system of loops is 
scale invariant, the above scaling behavior can be inferred 
from a dimensional analysis [7]. 

Now, the scale invariance of the above number den- 
sities, together with the self similar scaling of the loop 
length, can be used to obtain the scaling assumption 
n{£) ^ l^'^ for the number density of loops with length 
£, see Refs. [7I18I9J . defining the exponent r = 1 -I- djdf. 
To illustrate this it is useful to consider the number of 
loops with size from R^ to Rs -I- dRs , given by dN — 
n{Rs) dRs- Using the leading order of the relation be- 
tween the loop length and loop size, i.e. £^Rs and hence 
dRs^ £~^~^^^'^-^ d£, one can transform dN by a change of 
variables Rs -^ £ to yield 

dN - Rj^'^+^'> dRs - £-(1+''/^/) d£ - n{£) d£. 
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Fig. 6. Results obtained for the pdf Pl{£) of the con- 
figurational energy per lattice site e. (a) Data obtained 
for DPLs on the 2d square lattice, (b) data for DPLs on 
the 3d simple cubic lattice. In either case, the shape of 
Pl (e) compares well to a Gaussian distribution (indicated 
by the dashed line). The insets illustrate the average and 
the instance to instance fiuctuation of e as function of the 
system size L. 



The scaling assumption for n{£) above fits the data 
quite well, where, for the setups considered here, we yield 
the estimates of r listed in Tab. [TJ The data for the 2d-sq 
and 3d-sc DPLs is shown in Figs. [SIa),(b), respectively. 
Finally, note that the values of t and df listed in Tab. [1] 
where obtained independently and agree with the scaling 
relation above within errorbars. 



O. Melchert, A. K. Hartmann: Configurational statistics of DPLs and FPLs in the NWP model 



Note that, for the Sd-sc DPL setup, if one considers 
all loops (not only the non-spanning loops) in order to 
compute the distribution n{i), one finds a crossover to the 
scaling behavior n{£)^£~^ in the limit of long loops (i.e. 
£> 10^ for the 3d-sc lattice with L = 48, see Fig.Efb)). This 
behavior is a direct consequence of the periodic boundary 
conditions that allow for long loops that wind around the 
lattice, i.e. those loops with a spanning index j>0. Such 
a crossover behavior was also observed within simulations 
of other string bearing models [33134] , and, by analytic 
means, for a simple model of string formation on periodic 
lattices [T^ . 



3.2 Energetic properties 

Scaling of the configurational energy: The disorder av- 
eraged total weight {£) of the loop configuration is an 
extensive quantity for which we found a clear asymptotic 
scaling that allows us to estimate numerical values for 
(e) = {\£\/L'^) as listed in Tab. [T] Here, we found no clear- 
cut deviations from the scaling behavior that might allow 
for an analysis of corrections to scaling. In this regard, 
for the smallest 2d-sq DPL system considered we found 
e = 0.5965 (L = 12), while for the largest system we ob- 
tained £ = 0.5958 (L = 256). For the 3d-sc DPL setup we 
similarly find e = 0.84447 (L = 10) and e = 0.84446 (L = 48) . 
Moreover, the pdf Pl{^) compares well to a Gaussian 
distribution, as shown for 2d-sq and 3d-sc DPLs in Figs. 
[nKa),(b), respectively. We further analyzed the instance to 
instance fluctuations {6e) — [(e^) — (e)^]^/^ (inset of Figs. 
[6i;a),(b)), for which we estimate (5e) =0.680(4)L-0-999(2) 
(2(i-sq DPLs) and (Se) ^Om{l)L-^-^°*^^^ (3d-sc DPLs). 
We therefore deduce the asymptotic scaling (Se) ^L^"^/"^. 
Note that in terms of the extensive configurational energy 
£", the latter reads {5£) ^ L'^^'^ and is thus in agreement 
with the result for finite dimensional spin glasses |35j . As 
evident from Fig. [HI the data sets for different system sizes 
show a nice data collapse if rescaled according to Pl (e) = 
i-rf/2p^[(e_(e))/irf/2]. 

Considering DPLs for 2d square and honeycomb lattice 
graphs, one might speculate whether the respective esti- 
mates of (e) reduce to a unique value, apart from purely 
geometric factors. In this regard, assuming that e depends 
on the fraction of occupied lattice sites pcov as well as the 
coordination number z that characterizes the particular 
lattice geometry (2d-sq: z — A; 2d-hy: z — 3), we find that 
the simple expression (e) =0.182 -z-pcov fits the respective 
values of (e) (listed in Tab. [1]) quite well. 



Scaling properties of the loop weight: Further, we in- 
vestigated the scaling of the average loop weight {ojc) 
as function of the loop spanning length R by consider- 
ing the full ensemble of loops collected at the largest sys- 
tem size simulated. In either case we find a good fit to 
the scaling form (ojc) ~ i?"(l -I- cR~^), where we yield 
the exponents a = 1.759(7) (2.05(7)) and /3 = 0.84(4) 
(0.58(7)) for 2d-sq (3d-sc) DPLs (results for the 3d-sc 
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Fig. 7. Scaling of the probability density function (pdf) 
Pr{ujc) for the loop weight cj£ of loops with spanning 
length R. (a) the main plot shows the scaling of the pdf for 
2d DPLs and FPLs on a semi-logarithmic scale, while the 
inset shows the same data on a linear scale, (b) scaling be- 
havior of the pdf for 3d DPLs, where the scaling function 
P2[\^c\/{^c)r\ (solid line) compares well to a log-normal 
distribution function as discussed in the text. The inset 
shows the scaling of the average loop weight as function 
of the spanning length R. 



data are shown in the inset of Fig. [7l[b)). For the 2d-sq 
FPL setup and considering a similar scaling expression, 
we yield a =1.75(2) and /3 = 0.87(17) (here, a more basic 
scaling function yields (ljc) ^ {R + 1.3(2))i-^4(i)). Note 
that a scaling according to the expression above, together 
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with the numerical values of the respective exponents, in- 
dicates that {ijjc)^t- Moreover, the root mean square de- 
viation {SbJcj = [{ijj'c) — {^jJcj'^Y^'^ exhibits the same asymp- 
totic scaling as the average loop weight, i.e. {Sujc) ^ R°' 
for sufficiently large values of R. However, for small val- 
ues of R there are also deviations from a pure power law 
scaling that might be accounted for by using a scaling ex- 
pression of the form {6loc) ^ {R + AR)" . E.g., for 2d-sq 
DPLs (FPLs) we find Zii? = 4(1) (12.5(3)) and a' = 1.75(2) 
(1.77(2)). 

In either case, the pdf Pr{ujc) of the loop weight for 
loops of different spanning length R exhibits the scal- 
ing behavior Pr{ll>c)'^ {SuJc)B^P2[i(^c - {^c)r)/{Sujc)r]- 
Comparing 2d-sq DPLs and FPLs, we find that both se- 
tups appear to have the same scaling function P2[-], i.e. 
after rescaling both data sets displayed in Fig. [Tlja) fall 
onto a unique mastercurve. 

For the 3d-sc case, we observe that {cjJc}r ^ {Sujc}r 
already for small values of R. Thus, the data fits 
a slightly simpler scaling assumption, i.e. Pr{i-^c) = 
{ujc)^^P2[\^c\/ {^c)r\^ where the scaling function P2[x\ 
can further be parameterized by the log-normal distribu- 
tion function P2[x\ = exp[— (log(a:) — ^)^/(2(T^)]/(xo-\/27r) 
with parameters n = —0.089(1) and a = 0.211(1), illus- 
trated by the solid line in Fig. [Tfb). 



4 Conclusions 

In the presented study, we performed simulations to char- 
acterize geometric and energetic properties of densely and 
fully packed configurations of loops in the negative- weight 
percolation model. Therefore, we considered a Gaussian 
distribution of the edge weights and 2d square/honeycomb 
and Bd simple cubic lattice graphs. Further, so as to obtain 
configurations of minimum weight loops, we used a map- 
ping of the NWP model to a combinatorial optimization 
problem that can be solved by means of exact algorithms. 
In subsequent FSS analyses we observed that the scaling 
behavior of observables is best represented by considering 
scaling expressions that also account for small power law 
corrections. 

Regarding the 2d setups, we found that the configura- 
tions of loops for both, DPLs and FPLs, are characterized 
by the same values of the geometric scaling exponents df 
and r. Moreover, the numerical value of df found here co- 
incides with the one measured for the smart kinetic walk 
(SKW) [33]. A loop-forming version of the SKW in two 
dimensions traces out the external perimeter of critical 
percolation clusters and was also used to model ring poly- 
mers at the ^'-point (which presumably is in the same 
universality class as the 0-point) on the honeycomb lat- 
tice [37] ■ Also, it compares well to experimental results on 
the scaling of the radius of gyration of 2d polymeric chains 
[55] . In that study, polymethylmetacrylate (PMMA) was 
spread onto an air-water interface at 16.5°C that acted 
as a 0-solvent for the respective polymer material. Fur- 
ther, the exponents df and r agree with those measured 
for fully packed loops on two discrete interface models 



termed "random manifold" and "random elastic medium" 
|23| . Finally, upon analysis of the pdfs that describe the 
scaling of the loop length as well as the loop weight, we 
found a universal behavior regarding the different 2d se- 
tups considered. 

For the Sd setup, we observed two phases of loops 
containing "small" non-spanning loops and "large" span- 
ning loops, respectively. Both phases are characterized by 
different geometric properties. Such a phenomenon was 
also observed for other string-bearing models, as e.g. cos- 
mic string networks [33 34 or more general models of 
string formation on periodic lattices [19118] . Regarding the 
small loop phase we conclude with a fractal dimension 
d/ sa 2.01(4). As discussed earlier, this scaling exponent 
indicates that the upper critical dimension for the DPLs 
takes the value d^^^ — Z. In contrast to this, earlier stud- 
ies revealed that the upper critical dimension of the NWP 
model at the critical point where percolating loops first 
appear in the limit of large system sizes is du = 6 [22] . As 
for usual self- avoiding lattice curves for d>d^^^, one can 
further expect to find d/ = 2. This we confirmed for DPLs 
on 4d hypercubic lattice graphs, where we obtained the 
estimate df = 2.00(5). Further, the large loop phase is a 
direct consequence of the periodic boundary conditions of 
the simulated systems. Upon analysis, we found that ap- 
proximately 90% of the sites on the 3d simple cubic lattice 
are visited by loops and large loops comprise about 91% 
of the loop segments on the lattice. Hence, these loops fill 
extensive parts of the "simulation box" and are strongly 
affected by the finite size of the accessible volume. Within 
the large loop phase we found evidence that the scaling 
properties of the loops further depend on their precise 
topology. 

Finally, note that the algorithmic procedure presented 
in section [2] is not limited to a particular disorder distri- 
bution or lattice setup and thus allows for the analysis of 
configurations of DPLs and FPLs for a multitude of lattice 
geometries and lattice dimensions. 

OM acknowledges financial support from the VolkswagenS- 
tiftung (Germany) within the program "Nachwuchsgruppen 
an Universitaten" . The simulations were performed at the 
GOLEM I cluster for scientific computing at the University 
of Oldenburg (Germany). 
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